Equilibrium Properties of a Trapped Dipolar Fermion at Finite Temperatures 
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We study the equilibrium properties of a dipolar Fermi gas at finite temperatures. We introduce 
a variational ansatz for the phase-space distribution function that can describe the deformation in 
both real and momentum space. The effect of dipole-dipole interactions on the thermal equilibrium 
is discussed with particular emphasis on the deformation in momentum space. We examine the 
stability of the system by varying the temperature, trap aspect ratio, and the dipole moment. In 
addition, we discuss how the deformation in both real and momentum space can be observed in the 
high-temperature regime. 



Interest in dipolar gases has been growing since Bose-Einstein condensates (BECs) of 52 Cr atoms, which have 
large magnetic dipole moments, were experimentally observed The anisotropic and long-range nature of the 

dipole-dipole interaction confers rich properties to both the equilibria and dynamics of dipolar gases. There have 
been a number of theoretical studies of dipolar BECs, which have investigated their ground states [1, HI, collective 
oscillations [H, 0] , and their properties in optical lattice potentials [1, Q . The ground state [1(1 [HI and stability 



12 1 of dipolar Fermi gases have been studied. In addition, there have been theoretical studies of expansion 



14[, collective oscillations [15|, and the supcrfluid phase [lg-|l9|. On the other hand, several experiments proceed 



energetically toward the realization of creating heteronuclear polar molecules whose large electric dipole moments 
result in strong dipole-dipole interactions. Rydberg atoms have also been attracting attention due to their large 
electric dipole moments. 

Most studies of dipolar Fermi gases have focused on the zero-temperature properties of these gases and there have 
been few theoretical studies on finite-temperature gases [l9rl2ll |. Despite many groups conducting experiments, no 



groups have succeeded in cooling polar molecules to the quantum-degenerate regime [221 ] . It is thus important to 
investigate the temperature range in which the dipole-dipole interaction has appreciable effects. For this reason, we 
concentrate on the properties of dipolar Fermi gases at finite temperatures. 

Unlike dipolar BECs, dipolar Fermi gases do not interact via s-wave collisions and they have both Hartree direct and 
Fock exchange energies of the dipole-dipole interaction in the mean-field description, which reflects the antisymmetric 
many-body wave function. Miyakawa et al. introduced a variational Wigner function to examine the ground state of 
the system at zero temperature and they demonstrated that the anisotropic nature of the interaction causes Fermi 



surface deformation through the Fock exchange energy [llj . In the present paper, we generalize this variational method 
to finite temperatures and discuss the effect of temperature on the equilibrium properties. In particular, we focus on 
the deformation of the distributions in real and momentum space and investigate the temperature dependence of the 
instability of dipolar Fermi gases against collapse. 

We consider trapped dipolar fcrmions. The dipoles are assumed to be polarized along the z axis due to an external 
electric field. In the second quantized form, the Hamiltonian for this system is given by 



H = J dr¥ (r) 



■£v 2 + Wr) 



*(r) 



+ \J dv j dv'¥ (r) ¥ (r') V dd (r - r') <F (r') <F (r) , (1) 

where ^ (r) is the Fermi field operator, and the hat indicates a second quantized operator. The first term of Eq. 
(P) describes the Hamiltonian of a single particle in a harmonic trap V tra p (r) = (m/2) [w 2 (x 2 + y 2 ) + u 2 z 2 ~\ , where 
m is the particle mass. The second term describes the two-body interaction Hamiltonian of the dipole-dipole force, 
where Vdd (r) = (d 2 /r 3 ) (l — 3z 2 /r 2 ). d is the coefficient of the dipole-dipole interaction and d 2 = p 2 /Atteq, where p 
is magnitude of the electric dipole moment [ill . [l5| and eo is the electric permittivity of a vacuum. We introduce the 
Wigner distribution function: 

W(p,r) = J dr'e- l v r ' /h {¥{v-Y'/2)i>{v + v'/2)). (2) 

The density distributions in real and momentum space in terms of W (p, r) are given by n (p) = J dr W (p, r) and 
n (r) = (2irh)~ 3 J dp W (p, r), respectively. 
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We consider the thermal equilibrium of dipolar fcrmions trapped in a harmonic potential. The system is not globally 
stable when a dipole-dipole interaction is present because the interaction is partially attractive and causes collapse of 
the gas. However, a metastable state exists at finite temperatures under certain conditions. To find this metastable 
state, we introduce a variational Wigner distribution function that is analogous to the one introduced in Ref. [llj ]. 
We assume the Maxwell-Boltzmann regime at relatively high temperatures, namely 



W (p, r) = exp 




Mo 



/k B T 



(3) 



1 /3 

where p 2 = x 2 + y 2 , p 2 = p 2 +Py, and oj = (cj 2 cj z ) . Here, the positive parameters a and j3 represent deformations 
of density distributions in momentum and real space respectively, and A describes compression of the dipolar gas, as 
defined in Ref. [ll|. In addition to these parameters, we introduce a new variational parameter 9 that characterizes 
compression in momentum space. In addition, the chemical potential is determined by the number constraint N = 
J dr n (r) , which gives 



From the normalization condition 



N 



dr 



dp 
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-W(p,r), 
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the chemical potential po is determined by Eq. Q. Under this ansatz, the density distribution in momentum space 
is given by 



n (p) = N 



2/12 X 3/ 2 



mk B T 



exp 
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(6) 



Thus, the aspect ratio of the momentum space distribution, which is the ratio of the root-mean-squarc momentum in 
the p z direction to that in a direction in the p x ~Py plane (we choose the p x direction in the following) yl (p 2 ) / (p 2 ) , 
becomes a~ 3 / 2 . Similarly, we obtain the density distribution in real space as 



V kBT 
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(7) 



leading to the aspect ratio of the real-space distribution (i.e., the ratio of the root- mean-square radius in the z direction 
to that in the x direction) being a/ (z 2 )/(x 2 ) = (3 3 > 2 . We note that, from Eq. (jTj) , A > 1 (A < 1) corresponds to 
compression (expansion) of the gas in real space, which is a consequence of the effective attraction (repulsion) of the 
dipolc-dipolc interaction. In addition, comparison of Eq. ([6]) with Eq. (JT)) reveals that 9 plays the same role in real 
space as A plays in momentum space. 

To find a metastable state at a temperature T, we look for a local minimum of the Helmholtz free energy 



F = E -TS. 



(8) 



Here, the total energy E is the sum of the kinetic energy Ex, the trapping potential energy Ey, the Hartree direct 
energy Eh, and the Fock exchange energy E ex : 



E = E K + E v + E H + E ex , 
where the four contributions are given in terms of the Wigner distribution function as [Til ] 
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(10) 

(11) 
(12) 

(13) 



3 



In Eq. (IT3l) . we have introduced the center of mass coordinate R = (r + r')/2 and the relative coordinate 
s = r — r'. In deriving Eqs. ([12"]) and ([i"3]) . we used the mean-field decoupling (r) ^ (r') ^ (r') (r)) ~ 
(# t (r) * (r))(* t (r') * (r')) - (r) ^ (r'))(*t (r') * (r)}. We note that the minus sign in Eq. ([13]) arises from 
the Fermi statistics. The entropy S is given by 

S = -k B fdr f -^{W(p,r)\n[W(p,r)] + [l-W(p,r)}\n[l-W(p,r)}} 
J J (27m) 

~ -k B fdr [ (p, r) In [W (p, r)] . (14) 

J J (27m) 

In the above expression for the entropy, we have assumed W <C 1, which is consistent with the Maxwell-Boltzmann 
distribution. 

Using the variational function (j3]) in Eqs. (fTU|) - (fl"4"]) . we can express the free energy in terms of the variational 
parameters: 

{1 1 C\ 3 



-3 + In 



N(\9ujh) 3 - 3\n[k B T]^Nk B T, (15) 



where C = (d 2 rn z / 2 us' i iV) / (3 ■ 2 2 7r 1 / 2 ). In the above equation, we have defined the scaling functions / (a) = 2a + l/a 2 

and g (f3) = 2(3 Q //3 + /3 2 //3 2 , where (3 = (oj p /oj z ) 3/2 . The deformation function / (a) is defined in Refs. [TJ d- The 
first and second terms in the first line of Eq. (| 15[) are the kinetic energy and the potential energy, respectively. We 
note that the kinetic energy is inversely proportional to the square of 9. Thus, an increase in 9 indicates a reduction in 
the effective temperature T* = T/9 2 associated with the kinetic energy. The third and fourth terms are the Hartree 
direct energy and the Fock exchange energy by the dipolc-dipolc interaction, respectively. Note that the deformation 
parameters a and /? contribute separately to the internal energy, E. The momentum space deformation parameter a 
appears only in the kinetic and Fock exchange energies, whereas the real space deformation parameter (3 appears only 
in the potential and Hartree direct energy. The last three terms in Eq. ([15]) express the entropy, which is independent 
of a and /3. 

The equilibrium solution can be found by minimizing the free energy F with respect to the four variational param- 
eters, a, (3, A, and 9. To be a local minimum of F representing a metastable state, the solution must be a convex 
downward point in the four-dimensional space of the variational parameters. This requires the following conditions: 

< a < 1, ri „x 

/?o < A (16) 

From the conditions (|16[) , wc find that the momentum space distribution is always elongated in the p z direction The 
real-space distribution tends to be stretched in the z direction and the aspect ratio of the real-space distribution 
is larger than that for the non-interacting case. These results arc consistent with those for the zero-temperature 
case fill ]. 

We find the equilibrium solution by minimizing Eq. (|15[) numerically by varying three system parameters (temper- 
ature T, trap aspect ratio /3o, and the magnitude of the electric dipole moment p) for m = 100 a.u.m., u = 2tt x 10 2 
Hz, and N = 10 4 . 

Figure [T] shows a plot of the variational parameters as a function of temperature in units of ideal gas Fermi 
temperature Tg = (6A^) 1/3 huj ~ 1.88 nK for /3 a = 0.5 (solid line), /3 = 1.0 (dashed line), and /3 Q = 2.0 (dash-dotted 
line) and for when the electric dipole moment is fixed at p = 1 Debyc. Figure [TJa) indicates that the momentum 
distribution becomes more elongated in the dipolar direction as the temperature decreases. The deformation in the 
momentum distribution arises from competition between the kinetic energy, which favors an isotropic momentum 
distribution, and the exchange energy, which favors an anisotropic momentum distribution. As the temperature 
decreases, the kinetic energy decreases while the interaction energy increases, leading to a large aspect ratio in the 
momentum distribution. Figure Eh) reveals that a lower temperature gives rise to a larger deviation of the aspect 
ratio in real-space from the non-interacting case. It also shows that the shift of /3 is largest in the isotropic trap. 
FigureQTc) shows that the magnitude of the deviation from 1 in A is larger at lower temperatures. In addition, we have 
A > 1 for /3 > 1 (cigar-shaped trap) and A < 1 for (3 < 1 (oblate trap) . The Hartree direct energy becomes attractive 
(repulsive) in the former (latter) case due to deformation of the spatial density distribution mainly determined by 
the trap geometry, which leads to compression (expansion) of the gas cloud for f3o > 1 (/?o < 1). Finally, Fig. Old) 
reveals that the dipolar Fermi gas is compressed in momentum space irrespective of the trap aspect ratio because only 
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FIG. 1: Variational parameters as functions of the temperature T for a fixed dipole moment p — 1 Debye for a trap asymmetry 
j3o — 0.5 (solid line), /3o = 1.0 (dashed line), and /3q — 2.0 (dash-dotted line). 



the kinetic energy and entropy depend on 9 (Eq. (| 15[) ) . In addition, the momentum distribution is not compressed 
when there is no exchange energy. Whereas the interaction effect increases at lower temperatures in all cases, the 
deformation in both real and momentum spaces are invisibly small for T > 2Tp and p = 1 Debye. 

Figure [5] shows the variational parameters as a function of temperature for p = 1.0 Debye (solid line), p = 3.0 Debye 
(dashed line), and p = 5.0 Debye (dash-dotted line) for an isotropic trap /3o = 1- In all cases, the deviation of the 
variational parameters from those of the non-interacting case becomes more pronounced at lower temperatures and at 
higher electric dipole moments. For large electric dipole moments, we plotted the optimized values for T > T c , where 
T c is the critical temperature below which the variational free energy (|15p has no local minimum and the system 
becomes unstable to collapse. The critical temperature for this collapse for p = 5.0 Debye is T c — 2.2hT F . On the 
other hand, systems with p = 1.0 Debye and p = 3.0 Debye are stable in the temperature regime T > 1.5Tp. These 
results reveal that the system is unstable even in the high-temperature regime for large electric dipole moments. 
Figure [2ja) and (b) shows that the aspect ratios in both momentum and real space are sufficiently large to observe 
deformation effects in experiments for polar molecules with large electric dipole moments. 

Figure [3] shows the stability diagram as functions of the electric dipole moment p and trap aspect ratio /?o for 
different temperatures. This result indicates that the critical dipole moment p c increases drastically as the trap 
becomes more oblate. This is consistent with the result for the zero-temperature case [ll|. The unstable region 
expands with decreasing temperature. From Fig. [21 we presume that the critical temperature for collapse instability 
increases as the electric dipole moment increases. This behavior is confirmed by Fig. SI which shows the critical 
temperature T c for /3o = 0.5, /3o = 1-0, and /3q = 2.0. Figure [S] shows the critical temperature as a function of the trap 
aspect ratio for p = 1.0 Debye and p — 2.0 Debye. This result indicates that the critical temperature T c increases 
more as /?o increases. When p= 1.0 Debye, the system is always stable within the parameter range shown in the 
figure (Po < 40). Figure \S\ suggests that we should be careful not to enter the unstable region when the temperature 
is reduced. Even if the system is initially in a stable region at a high temperature, the system may (depending on 
the trap aspect ratio) become unstable on cooling before reaching the quantum degenerate regime. In addition, we 
see that the local minimum always disappears when A exceeds a critical value. This suggests that the instability first 
occurs in real space and a dipolar Fermi gas collapses when A reaches the critical value. 

Finally, we discuss the stability properties in connection with experiments. The JILA experiment 22j revealed that 
the mass of the polar molecule is about m— 128 a.m.u. and that the electric dipole moment is about p= 0.566 Debye. 
Taking these values along with u> = 2tt x 10 2 and N = 10 4 , it is always possible to find a stable dipolar Fermi gas for 
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FIG. 2: Variational parameters as functions of the temperature for a fixed trap asymmetry /3q = 1.0 for a dipole moment 
p = 1.0 Debye (solid line), p = 3.0 Debye (dashed line), and p = 5.0 Debye (dash-dotted line). 
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FIG. 3: Critical value of the dipole moment p c as a function of the trap aspect ratio for T/T F = 1.5 (open circles), T/T F = 2.0 
(filled circles), and T/T F = 3.0 (crosses). 




FIG. 4: Critical temperature as a function of the dipole moment for trap aspect ratios /3o = 0.5 (open circles), /3q — 1.0 (filled 
circles), and /3o = 2.0 (crosses). 
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FIG. 5: Critical temperature as a function of the trap aspect ratio for dipole moments p = 3.0 (filled circles) and p = 5.0 
(crosses). 



any temperature T > T F in the trap aspect ratio regime fio < 40. Since, according to the results in Ref. [TTj], a g as 
of polar molecules at zero temperature with the same system parameters becomes unstable, the critical temperature 
should lie in < T < T F . However, the dipolar Fermi gas becomes unstable at extremely large trap aspect ratios. 
For ui = 2n x 10 2 , N = 10 4 , p = 0.566 Dcbye, and m — 128 a.m.u., the critical trap aspect ratios for T = 1.5T F and 
T = 2.0T F are /3g = 62.0 and /3g = 72.5, respectively. 

We discuss how the deformation effects of momentum-space and real-space distributions can be observed in the 
polar molecules used in the JILA experiment j22l |. For the same parameters given in the previous paragraph except 
for the frequency of the harmonic oscillator potential u, increasing u causes the dipole-dipole interaction effects to 
become very pronounced. This is because the density of the gas trapped in a tight potential becomes high and 
the interaction becomes more effective. For tu — 5 x 2ir x 10 2 Hz, /3o= 1.0, and T = l.bTp, we have aspect ratios 
\/(Pz)/ (Px ) = a -3 / 2 ~ 1.130 in the momentum distribution and ■J (z 2 ) / (a; 2 ) = /3 3 / 2 ~ 1.127 in the spatial distribution. 
For u = 6 x 2?r x 10 2 Hz, j3 = 1.0 and T = 2.0T F , we have a~ 3/2 ~ 1.106 and /3 3 / 2 ~ 1.103. Thus, the deformation 
effects arising from the anisotropic nature of the dipole-dipole interaction are detectable in the present experimental 
conditions. 

In conclusion, we have used a variational method to study the equilibrium properties of a dipolar Fermi gas at finite 
temperatures. As at zero temperature, the anisotropic nature of the dipole-dipole interaction leads to deformations 
in momentum and real space and the partial attraction of the interaction causes instability of the gas to collapse. In 
addition, we found that the dipolar Fermi gas becomes compressed in momentum space as the electric dipole moment 
increases. We found that the deformations in both momentum and real space can be observed in the high-temperature 
regime with a large electric dipole moment and a high trap frequency. In addition, we found that the stable region 
expands at finite temperatures. These results will be useful when cooling polar molecules in experiments. We hope 
this study will stimulate further experiments on dipolar Fermi gas at finite temperatures. 

In future studies, we intend to study dynamics of a dipolar Fermi gas at finite temperatures, such as the expansion 
dynamics and collective oscillations. In particular, expansion dynamics is important since it can provide direct 
information of the deformation of momentum distribution. Sogo et al. have studied the expansion dynamics at zero 
temperature and found that the dipolar Fermi gas expands along the dipole moment direction irrespective of the 
trap aspect ratio fl5j . We will study the temperature dependence of the expansion dynamics in the presence of the 
dipole-dipole interaction. 

Note added. Recently, we became aware of two works of equilibrium properties of dipolar Fermi gas at finite 
temperatures [2fj|. [2l|. 
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